Asymmetric Nuclear Matter in the Relativistic Brueckner 

Hartree Fock approach 

F. de Jong and H. Lenske 
Institut fiir Theoretische Physik, Universitdt Giessen, 35392 Giessen, Germany 

(February 9, 2008) 



Abstract 

We present a calculation of asymmetric nuclear matter properties in a rela- 
tivistic Brueckner Hartree Fock framework. Following other calculations the 
components of the self-energies are extracted by projecting on Lorentz invari- 
ant amplitudes. It is shown that for asymmetric nuclear matter one needs a 
sixth invariant. We present a set of invariants which in the limit of symmetric 
nuclear matter reduces to the conventional set. We argue that the existence 
of such a properly behaving set is also crucial for the application of the pro- 
jection method in symmetric nuclear matter. Results for the equation of state 
and other observables are presented. Special attention is payed to an analysis 
in terms of mean-field effective coupling constants. Apart from the usual ones 
we also find significant strength in the isovector scalar channel, which can be 
interpreted as an effective (5-meson. 
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I. INTRODUCTION 



An important step in understanding the binding properties of nuclear matter was made 
with relativistic approaches to nuclear interactions. They were capable to resolve to a large 
extent the longstanding problem of the underestimation of the binding energy of nuclear 
matter found in non-relativistic calculations. Conceptually, the relativistic theory provides 
a very natural link to modern boson exchange interactions and - at least in principle - relates 
in-medium interactions to free space nucleon-nucleon (NN) scattering. 

Applications have mostly been concentrating on the description of symmetric nuclear 
matter Although being different in theoretical and numerical details the various 

approaches are in rather close agreement on the bulk properties of nuclear matter. Inves- 
tigations of asymmetric nuclear matter, however, are rare, most of them of a very recent 
date One reason is certainly that up to very recent years the data base on asymmetric 

systems was rather limited because stable nuclei range around Z/A ~ 0.5 with a rather small 
bandwidth of variations. Asymmetric matter was mainly of theoretical interest, except for 
astrophysical studies, as the structure of neutron stars The situation has changed com- 
pletely with the new radioactive beam facilities producing now new isotopes in a sufficiently 
large amount for detailed studies. The isospin degree of freedom in the nuclear equation 
of state can be investigated for the first time over a much larger range of proton-neutron 
asymmetries as before. At present, the largest asymmetries produced experimentally are 
obtained for light dripline nuclei {Z/A = 0.2 ■ ■ -0.3 for the heavy isotopes of He to oxygen 
nuclides) which is about a factor of 2 less of the asymmetries found in stable nuclei. In 
the near future an increasing amount of data for nuclides all over the mass table can be 
expected. 

It is obvious that the properties of asymmetric matter are of high interest for nuclear 
structure at the proton and neutron driplines and vice versa. The isospin dependence of 
the nuclear equation of state is to a large extent unexplored. Extrapolations of the nuclear 
mass formula into the regions of new isotopes persistently result in strong deviations from 
measured binding energies. The same is found for structure calculations with seemingly well 
established effective interaction models as the non-relativistic Skyrme energy functional. 
Although the parameters of Skyrme interactions are derived empirically and the various 
model interactions reproduce the properties of stable nuclei very well the predictions for 
nuclei far off stability are found to differ drastically. Similar observation are also made for 
relativistic mean-field theories. Thus, the question arises to what extent the asymmetry, i.e. 
isospin dependent contributions to the Bethe-Weiszacker mass formula are really understood. 

It is apparent that systematic studies of asymmetric matter are required. In this work 
we describe asymmetric nuclear matter in the relativistic Brueckner-Hartree-Fock (RBHF) 
approach. Compared to symmetric matter the theoretical and numerical efforts, respectively, 
are much larger because protons and neutrons are occupying different Fermi spheres. This 
implies different effective masses for protons and neutrons and the two sorts of particles cover 
a different range of off-shell momenta. Hence, it can be expected that beside the isospin 
also the momentum dependence of self-energies will be of importance. As noted before, 
asymmetric matter requires to solve a system of coupled equations for the set 7Vi,r2,r3,T4 of in- 
medium proton-neutron T-matrices, where r = ibl denote a proton or neutron, respectively. 
We perform our calculation in a plane-wave spinor basis in which we have neutrons and 
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protons as distinguishable particles. This implies that we also take into account partial- 
wave amplitudes that are zero in symmetric nuclear and neutron matter, where the masses 
are equal and isospin is a good quantum number. It is not clear whether other calculations 
take these amplitudes into account. Details of the theoretical formulation are discussed in 
Sect. 2. 

A central point of the discussion is the derivation of a new asymmetry dependent invari- 
ant. Traditionally, the positive energy-projected in-medium on-shell T-matrix is expressed 
in terms of five Lorentz invariant amplitudes (see Eq.(2.7)). This standard set of invariants 
is known to be not unique because each can be replaced by its derivative counterpart. A well 
known example is the ambiguity found by changing from the pseudo-scalar to the pseudo- 
vector vr — A^A^ coupling. A proper treatment avoiding such ambiguities is to solve the 
in-medium scattering problem for the full Dirac-space including also negative energy states. 
Such an approach was applied by Huber et al. p using the so-called A-approximation. In 
our approach we continue to use the projected T-matrix. To do so we need to introduce 
a sixth invariant since in asymmetric matter there are six independent helicity amplitudes. 
Clearly, the new invariant must be defined such that the description of symmetric matter 
is conserved. A promising guideline is to consider an amplitude which approaches zero for 
equal proton and neutron effective masses. Such an amplitude will vanish in symmetric mat- 
ter and thus fulfills the constraint mentioned above. In order to account for the differences 
in Tnp^np and Tnp^pn one actually has to define separate invariant amplitudes for the direct 
and the exchange channel. In the direct channel a vector-scalar and in the exchange channel 
a pseudovector-pseudoscalar combination is required. However, only the direct amplitude 
contributes to the proton/neutron self-energies. 

Results for the equation of state for nuclear matter at various degrees of asymmetry are 
discussed in Sect. 3. Particular emphasis is put on the momentum structure of self-energies. 
We define effective isoscalar and isovector self-energies which directly characterize the mean- 
field of asymmetric nuclear matter. The dependence of the equation of state on the isovector 
density is tested and compared to the prediction of the mass formula of being quadratic. 
The paper closes in Sect. 4 with a summary and conclusions. 

II. THEORETICAL FRAMEWORK 

The implementation of the RBHF model we use is identical to the one presented in Ref. 

an account of a calculation for asymmetric nuclear matter can be found in Ref. We 
will only briefly review the concepts of the relativistic Brueckner model, but will elaborate 
further on an omission in the calculation of Ref. [^]. 

The central quantity in the model is the self-energy, in the RBHF model it has a structure 
in the spinor representation, expressed by scalar and vector components: 

S(p) = S^(p)-7°S°(p)+7-pS". (2.1) 

Inserting this form of the self-energy into the Dirac equation, it is natural to define 'effective' 
quantities by 

m* = m^v + S^(p), pl=Po + TP{p), p* = {l + S^(p))p. (2.2) 
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This can be simplified further by dividing out the (small) H'" term 0,^: m* = 771^ + 
S*'(j9)/(1 + S^(p)). One thus finds a solution of the 'dressed' Dirac equation; the effective 
spinor: 



/ 1 

\Xr, (2.3) 
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with E* = yp^ + m*^. The self-energy is found to depend only weakly on the momentum 
at least inside the Fermi-sphere: in the practical calculation it is taken constant and set to 
the value at the Fermi-surface. This greatly simplifies the calculations. 

Obviously in asymmetric nuclear matter the Fermi-momenta of protons and neutrons 
are different, leading to different Pauli-blocking operators and corresponding neutron and 
the proton self-energies. This asymmetry is also reflected in the fact that one has three dif- 
ferent effective interactions, represented by three different T-matrices: Tpp^pp (proton-proton 
scattering), Tpn^pn/np (proton-neutron scattering) and Tnn,nn (neutron-neutron scattering). 
In the Brueckner scheme these are calculated by: 

Tnn,nn ^nn,nn j ^nn,nn9nnTnn,nn 

'-^pp,pp ^pp,pp ~l~ J ^pp,pp9ppTpp,pp 

Tnp,np ^np,np ~l~ j ^np,np9npTnp,np ~l~ j ^np,pn9pnTpn,np 

Tpn,np ^pn,np ~l~ J ^pn,np9npTnp,np ~l~ J ^pn,pn9pnTpn,np (2'4) 

Note the coupled channels for T^p, we also stress that Vnp,pn is not related by a Fierz- 
transform to Vnp,pn'- in we have neutron-proton vertices on which only the isospin-1 

mesons contribute. It immediately follows that the same holds for the T-matrices: Tnp,np is 
not related to Tnp^pn- The 2-body propagators 9ij are the usual ones as they appear in the 
three-dimensionally reduced Thompson equation: 



E* E*^-E*-E* + ie ^ ' ' 

where Qij is the angle averaged Pauli blocker of the intermediate momentum p, a function of 
the total momentum P and invariant mass \fs* , its parameters depend on the appropriate 
Fermi-momenta Pf,i/j- After obtaining a solution for the T-matrices, the self-energy is then 
calculated using a Brueckner-type expression: 
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Tr 



Tl9n + / Tr T^op , (2.6) 



where T^ stands for Tnp,np{0 = 0) — Tnp,pn{6 = tt). For the proton self-energy the expression 
is completely similar. We will specify these expressions further later on. For the moment we 
note that the integrations extend over the neutron and proton Fermi sea, respectively. This 
calculation is repeated until a self-consistent solution for the self-energy is found. With the 
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self-consistent self-energy one then can calculate quantities like binding energy and single- 
particle energies. 

In our implementation of the relativistic Brueckner model we calculate the self-energy 
by decomposing the various T-matrices in Lorentz invariant amplitudes . This boils down 
to transforming the spin-representation to spinor-representation. The contributions of the 
various amplitudes to the components of the self-energy are then easily found. In the case 
of identical particles (to be more precise, particles with equal mass), five amplitudes are 
needed for the decomposition of the on-shell T-matrix. We use a set of invariants including 
the pseudo-vector coupling. In the cm. system we then have: 

(p'|T(P, V^)|p) = ^r,(P, V^,J9',J9)<I>, 

i 

/, G{l,7^a'^^7V,/4}• (2-7) 

2m* 

With q the transferred momentum: q = p' — p. For the on-shell scattering matrix elements 
we have in the cm. frame \p'\ = \p\ and \/s* = E* + E*. To calculate the self-energy one 
needs the amplitudes for the direct {9 = 0,p' = p) and exchange {6 = T[,p' = —p) seperately. 
At these points the set of invariants is degenerate, this is solved by decomposing at small 
angles and extrapolating to 6 = 0, vr. In Ref. it is shown that this limit exists. 

It has often been noted this particular set of invariants is not unique. One can always 
replace an invariant with its derivative counterpart, we did this by using the pseudo- vector 
invariant instead of the pseudo-scalar invariant. In this specific case one finds the same co- 
efficient Fj since the T-matrix is on-shell, but the contribution to the self-energy is different. 
However, the contribution to the total energy (found by sandwiching the self-energy between 
effective spinors) is the same for a given T-matrix. In a self-consistent calculation there will 
be a difference due to the fact that the self-consistent self-energy will be different, and thus 
the T-matrix is different. One guideline to limit the freedom in the choice of invariants is to 
require that one uses the same vertices in the invariant amplitudes as appear in the OBE 
interaction. Since we have a pseudo-vector coupling for the pion, we thus also have to use 
the pseudo-vector invariant in the projection. 

To avoid this ambiguity one has to calculate the scattering in the complete Dirac space, 
as is done e.g. in Ref. [Q. For the unique decomposition one needs u*{p)T,{p)u*{p), 
u* {p)T,{p)v* (p) = v*{p)Ti{p)u*{p) and v* {p)T,{p)v* (p) . In our calculation we only have the 
amplitude between positive energy spinors (the mean-field generated by the self-energy). 
Two different decompositions of the self-energy S, S' need to give the same value for the 
mean-field: u* {p)T,{p)u* (p) = u* {p)T,' (p)u* (p) . One can easily check that up to 0{p^/m*'^) 
this implies for the components S'^ = — a, Sq = Eq — a and S^^ = — a/2m*. The 
'ambiguity' coefficient a will of course depend on density and momentum. Inserting this in 
the expression for the renormalized effective mass we see that this is shifted by a quantity 
much smaller than a. Typically, an error of 100 MeV leads to a shift of ~ 40 MeV in the 
renormalized effective mass. From the above it follows that any approximate method for 
the decomposition of the self-energy will still give fair results for quantities that involve only 
the mean-field, like the binding energy. 

In the case of particles having different masses, the helicity matrix elements cease to be 
symmetric under exchange of particles. Since the number of invariants needed to decompose 
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the T-matrix is equal to the number of independent hehcity matrix elements after exhausting 
all the symmetries, the loss of this symmetry implies that one needs additional invariants 
to decompose the T-matrix. Appearantly this was not considered in Ref. ||^. For the on- 
shell np T-matrix, needed for the calculation of the self-energy, six invariants are required. 
Although it is relatively easy to find additional invariants that span the six- dimensional 
space of the six independent helicity matrix elements, the problem is that we need to have 
a decomposition that reduces to the one used in the symmetric case. With this we mean 
that the coefficient of the additional, sixth, amplitude should approach zero in the limit of 
equal masses. Obviously this is not a trivial requirement. 

Since T„p „p and T^p^pn are different we can choose invariants for the exchange and direct 
amplitudes independently. We encountered a similar situation for the A^A scattering matrix 
|Tl|. For the present case we found the following invariants which satisfy the above criterium: 



dir 
'6 

^exch 



p')u;{-p) +[<(p')<(p)] u;{-p')^u;{-p) 



uimi'H^P)] k(-p')7'<(-p)] + [uliP'WK^P)] k(-p')7'^<(-p)] , (2.8) 



where we used standard definitions from scattering theory k = pi — which results in = 
El — E4, k = p + p'. The numerical indices are the particle numbers (1 and 2 are incoming, 
3 and 4 outgoing), the indices n,p denote a neutron and a proton spinor respectively. The 
amplitude in the exchange channel is antisymmetric under exchange of particles, even when 
the masses are equal (the other invariants are symmetric under exchange of particles). This 
does not mean that it is not a valid invariant, but it results in a vanishing coefficient for 
equal masses, which is exactly what we required. We were unable to find an invariant 
with this property for the direct channel, but the one given above does have a vanishing 
coefficient for equal masses. In both invariants the momentum k appears. This makes 
it impossible to interpret them as some sort of meson channel since a meson carries a 
momentum q = Pi — Pz- To us there is no straightforward interpretation of these invariants 
beyond that they are needed to project out the T-matrix in the np channel. This T-matrix 
spans up a six-dimensional space in spin-space, and we therefore need a six- dimensional 
basis in spinor space for the basis transformation from spin to spinor representation. We 
want to stress again that it is crucial for the description of symmetric nuclear matter that 
for asymmetric nuclear matter we have a set of six invariants that behaves properly in the 
limit of symmetric nuclear matter (equal proton and neutron masses). If we could not find 
such a set, the method of extracting self-energies in symmetric nuclear matter by means of 
invariant amplitudes would have an additional ambiguity since the self-energies would not 
be a continuous function of the asymmetry parameter. 

Inserting the decomposition in the expression for the self-energies we find: 



Sn(P/,n) 



p/,n d?q 1 



(27r)3 2EI 



E {fTr[{f + m:)f ]rj,,,,„„ - f (r + m:)/TLMn} 



+ (1^^ E [fTrU + ^l)rn^r,np - f {f + "^P/T^^,} . (2.9) 



Working out the Dirac algebra, we find the components of the self-energy expressed in terms 
of the various r"(s*, P). Taking the Fermi- momentum pf along the z-axis and transforming 
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to the c.m.-defining variables s* — (pi + and P = pi +p2 we have for the scalar and 
vector components of the self-energy at the Fermi-surface 



Jo Js* 



ds* 



* 2 

xfdV" -V -AT'" -12r* -LAV"--" I ~ Pf,n1f^ ypv X 

V dir,nn exch,nn exch,nn exch,nn ' exch,nn ' 2,771*^ exch,nnl 



n 



1 m* 

as 



Jo Js* 



^^^'Pf,n £^;,,„ + El^ 
\ dir,np ' * dir,np exch,np exch,np exch,np ' exch,np 



2(m; - p'l^.^q^) - m; + 



ds* 



min.n 



\ Hir rtrt p-rrh nrt f^nh nm f^nh rtri i—t.,. ^ ... O f^nh nrj / 



l!lAlL\ 

dir,nn exch,nn exch,nn exch,nn 2771* ^ exch,nn) 



n 



Jo Js* . 327r2 



771* 



xf— ir" 4-m* g'P p6 _|_ ps — ^r'' — 2r'''' 

V dir.np '"'p tti* dir.np ' exch.np exch.np exch.np 

2E* (mf - p%q.) - E*(mf - ) 



2 



JO Js 



ds' 



^^^'Pf,n e;^^^ + El 



* 2 M 

/^pi) -I- _ 2r^ — 2r"'' — ~ Pf,n^f^ T^pv \ I 

V dir,nn ' exch,nn exch,nn exch,nn 9™*2 exch,nn) ' 



Jo Js* . 327r2 



z 



xf— 4r^ 4tt7* ^-^'^ ~ ^^ r^ + —2r^ — 2r"" 

V dir,np "^p dir,np ' exch,np exch,np exch,np 

Hz 



'^PfA^p -P^f,n%)-<iz^ 



g is the momentum of the integrated particle in the nuclear matter rest-frame, found by 
applying the inverse transform to the set (s*, P) Also is the on-shell four-momentum of 
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the incoming particle, taken at the Fermi momentum, so p° ^ = E*. The integration hmits 
are given by s^^^^^, = {E;^.+E;^^^ J^-P^ with q^in^i = \P-Pf,i\ and s*^^,^^ = AE;J-P^. 
The index i stands for either a neutron or a proton; pj i are the respective Fermi-momenta 
and E*, = (p2 + mf)i/2. 

The coefficient ^lxch,np does not appear in the expressions for the self-energy: for 6^ = vr, 
k = and since the extrapolated T^^^f^ .^^ is finite (we checked this by renormalizing the 
invariant by dividing by it does not contribute to the self-energy. 



III. RESULTS 

In Fig. |I] we present the results for the equation of state for various values of the 
asymmetry parameter Z/A using the Groningen potential [Q. Following Ref. we used an 
'averaged' Fermi momentum defined by: 

P=^Pl = ^PL + ^P?,n. (3.1) 

In the calculation it turns out that the contribution of the sixth invariant is negligible, 
although it of course indirectly affects the results by influencing the values of the other 
coefficients. Overall the results look pretty much like the ones presented in Refs. [0-^]. In a 
way this is not surprising because the results for Z/A = 0.5 (nuclear matter) and Z/A = 
(neutron matter) are unaffected by the sixth invariant. The binding energy curves for 
intermediate values of Z/A will lie between these two extreme curves, which does not provide 
a large freedom. A more quantitative measure of the equation of state is the asymmetry- 
energy coefficient 04, multiplying the 4{Z — A)'^/A term in the semi-empirical mass-formula. 
As has been argued in one should use the binding energies at the various saturation 
points in the determination of this coefficient and we determine 04 by 

1 

«4 = lo^^biO k=Z/A,p=p,,iZ/A) ■ (3.2) 

Using this prescription we find for the Groningen potential 04 = 25 MeV. Recalculating it for 
the Bonn C potential [§] we find a sightly larger value of 28 MeV. In the literature one finds 
a rather wide range of phenomenological values. Older ones tend to give lower values, e.g. 
in |jl2| a value of 04 = 23.7 MeV is given, a more recent liquid drop model calculation cites a 



result in the range 27-30 MeV while a very recent work comes up with 04 = 32.65 MeV |]T3| . 
Our result also shows that extracting the asymmetry coefficient via fitting self-energies only 
for symmetric nuclear and neutron matter, as is done e.g. in Ref. gives much too large 
results. These authors find values of 35 MeV and higher. We also note that although for 
the binding energy we find a nice quadratic dependence on the asymmetry parameter, this 
is not the case for the other quantities like effective masses and self-energies (taken at the 
respective Fermi- moment a). There we find a dependence which is more linear in character. 

In studying the momentum dependence of the self-energies it is convenient to define 
isoscalar and isovector quantities. For a given observable O we set 

Os = liOn + Op), = \{0n- Op), (3.3) 
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In Fig. 0we show the isoscalar and isovector mean-field as a function of the relative momen- 
tum p/pfi calculated at the respective saturation densities for each asymmetry parameter. 
As in Refs. we define the mean-field as the shift in the pole of the nucleon propagator 
due to the medium correlations. We see that up to an asymmetry parameter of 0.3 the 
isovector mean-field is rather independent of Z/A. There are some minor differences, but 
the slope is essentially the same. Only at Z/A = 0.2 we find a deviation, this might be 
due to the fact that we calculate the mean-field at the saturation density, which is different 
in this case. The isovector mean-field shows more structure. We observe the intuitively 
expected increasing difference with increasing asymmetry (corresponding to a lower value of 
Z/A). However, there is no quadratic dependence on the density A significant dependence 
on the momentum is found for high momenta, while inside the Fermi-seaa smooth behaviour 
is obtained. 

We show a similar analysis for the isovector components of the scalar and vector self- 
energy. The isovector scalar and vector self-energies can be interpreted as a 'effective' 
(5-meson and p-meson strengths. Again we see a notable momentum dependence, similar to 
the one observed in the isovector mean-field. 

It is also interesting to analyze our results in terms of mean-field effective coupling 
constants. This is particulary useful in calculations of finite nuclei [^]. As in the analysis 
above, we have four channels, both a scalar and vector in either the isoscalar and isovector 
channel. For this analysis we ignore the momentum dependence and calculate the self-energy 
of both the neutron and proton at the average Fermi-momentum as defined in Eq. p.l| . The 
effective coupling constants are then defined by: 



9iy ^ 


1 Kipf) + mpf) 




2 Pn + Pp 




1 + J^liPf) 




2 Pi + p; 




1 Kipf) - mpf) 




2 pi-pp 


= 


1 KiPf) - mpf) 




2 Pi -Pi 



(3.4) 

where pi and are the usual scalar and vector densities from relativistic mean-field theory 
]TP| . The results are presented in Fig. ^. As one intuitively expects, the scalar and vector 
isoscalar channels carry the largest strength. The dependence on the asymmetry parameter 
is small but we do observe a strong dependence on the density. For the densities presented in 
the figure, which range from half to double normal nuclear matter density, the dependence is 
predominantly linear in the Fermi momentum except for the isovector scalar channel. There, 
a significant dependence on higher orders mpf and the asymmetry is found. This is in agree- 
ment with the finding of Boersma and Malfiiet for their density dependent parametrization 
of the G-matrix. They strongly favour a density dependence of the coupling constants linear 
in pf as well [|T^, over other possibilities like the exponential dependence of Marcos et al. 
]r7| . The strength in the isovector channels is smaller although equal in both the scalar and 
vector channel. The significant scalar strength in the isovector channel suggests that for fits 



using only a limited number of meson, like the one of Ref. [17| one certainly needs to include 
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a 5-meson as well. For densities corresponding to pf > 0.24 GeV/c we again find negligible 
dependence on the asymmetry parameter and a linear dependence on the Fermi momentum. 
We find a positive sign for g'^'^, which is not the case when one uses the momentum depen- 
dence of the single-particle energy to extract the scalar and vector self-energy components 
p!8| . Note that since our decomposition of the self-energy is not unique the same holds for 
this analysis in terms of effective coupling constants. However, we think our approximate 
method still will be able to indicate trends. Still, both the neutron and proton self-energy 
will have an 'error' If these are very different due to e.g. a strong dependence on the 
momentum, one can obtain spurious results in the isovector channels as the negative sign of 
indicates when using the momentum dependence method. 

Below Pf ~ 0.24 GeV/c the density dependence changes and an enhancement of the 
isovector couplings is found at large asymmetry. Most probably this is related to the ap- 
pearance of bound states related to pairing, which typically appear at low densities in a 
Brueckner approach. In that sense the Brueckner scheme is an intermediate-density ap- 
proximation, losing its physical significance at low densities. At low densities the Brueckner 
independent-pair assumption ceases to be valid because of the onset of the pairing instability 
and the approach of the deuteron bound-state pole in the np-channel. For Z/A = 0.2 the 
proton density is already very low, and one starts to see the onset of pairing correlations. A 
similar behaviour can be observed in the density dependent parametrization of the former 



results of Ref. M by Marcos et al. |17 



The widely used relativistic mean-field theories (RMT) [|T9| , pO| are based on Lagrangians 
including o", uo and p mesons. Density independent meson-nucleon coupling constants ob- 
tained phemomenologically from fits to nuclear properties are used. In Tab. | our DB 
coupling constants, taken at the saturation density po=0.16 fm~^, are compared to standard 
density independent values from the a-u model (without scalar self-interactions) |T^. For 
comparison also the density dependent DB isoscalar scalar and vector coupling constants 
of ref. are shown which in []T3| were found to describe nuclear properties very well. In 
the isoscalar channels the two DB sets are in rather good agreement and both calculations 
reproduce almost perfectly the purely phenomenological RMT coupling constants. It is also 
worth noting that in the isoscalar channel the effective coupling constants at saturation den- 
sity are remarkably close to the 'bare' values as they appear in the interaction. This means, 
that around these densities the higher order effects generated by the Brueckner ladder are 
small. We therefore should find only small dependencies on e.g. the asymmetry, which we 
do observe. In the isovector channel the situation is totally different, there the effective 
coupling constants are much larger than the bare values. This implies that the (density 
and asymmetry dependent) higher order effects due to the Brueckner ladders are significant. 
Indeed, we see a pronounced dependence on asymmetry and density of the effective coupling 
constants of the effective isovector mesons. 

Most important and conclusive for field theoretical models are the results for the isovec- 
tor channels. The DB and RMT p coupling appear to be surprisingly close but is has to 
be noted that the isovector RMT coupling constants are not well determined. In fact, there 
exist other RMT parameter sets obtained from fits to finite nuclei with slightly different 
isoscalar couplings but about twice as large values for g^. From Tab. | it is seen that the 
isovector-scalar and the isovector-vector DB couplings are of comparable strength. This 
result gives clear evidence that the (5-meson should be included in effective field theories. 
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The scalar coupling of the S meson has important consequences for the dynamics in N 7^ 
Z systems. Because protons and neutrons will carry different effective masses they become 
dynamically distinguishable while the p meson leads primarily to energy shifts. Most impor- 
tant for applications to finite nuclei is that isovector effects in the spin-orbit potential will 
be enhanced as required by single particle spectra in charge asymmetric nuclei. Inclusion of 
the 6 meson might help to resolve the uncertainties in the RMT isovector sector. 



IV. CONCLUSIONS 

We presented a calculation of asymmetric nuclear matter in a relativistic Brueckner 
framework. To extract the components of the self-energies we used the method of projecting 
on Lorentz-invariant amplitudes. We found that for asymmetric nuclear matter a sixth 
invariant is needed on top of the conventional set of five invariants as for example used in 
Refs. [^0. This sixth invariant has to be chosen in such a way that for an asymmetry 
parameter close to 0.5 its coefficient has to vanish so that one recovers the usual description 
of symmetric nuclear matter in terms of five Lorentz-invariant amplitudes. We presented 
a choice of this sixth invariant that satisfies this criterium. With this model we performed 
calculations for several observables. The binding energy of the saturation point showed the 
expected quadratic dependence on the asymmetry parameter. For the asymmetry energy 
we found a value of around 25 MeV. We analysed our results in terms of mean-field scalar- 
vector isoscalar-isovector quantities. Apart from the expected scalar-vector strength in the 
isoscalar channel we also found significant scalar strength in the isovector channel. This can 
be interpreted as an effective 5-meson, which couples as strong as the effective isovector- 
vector p-meson. Finally, we found a density dependence linear in the Fermi momentum of 



the effective coupling constants, which was also observed in other calculations |T6 . 
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FIGURES 




FIG. 1. The binding energy as a function of the density for proton ratios ZjA ranging from 
(neutron matter) to 0.5 (nuclear matter). 
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p/p, (-) 



FIG. 2. The iso-scalar and iso-vector components of the mean-field for various asymmetry 
parameters. The mean-fields are evaluated at the saturation density for the respective asymmetry 
parameter. The momentum is given relative to the average Fermi-momentum as defined in the 
text. 



14 




FIG. 3. The isovector components of the scalar and vector components of the self-energy for 
various asymmetry parameters. These are evaluated at the saturation density for the respective 
asymmetry parameter. The momentum is given relative to the average Fermi-momentum as defined 
in the text. 
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FIG. 4. The effective coupling constants in the isoscalar scalar (ga), and vector (g^^) and the 
isovector scalar (gg) and vector (gp) as defined in the text for various asymmetry parameters and 
densities. 
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TABLES 



Meson 








Mass [MeV] 


a 


7.82 


7.51 


7.29 


550 


u 


12.10 


11.52 


10.84 


783 


P 


2.97 




2.93 


770 


5 


4.61 






983 



TABLE L DB coupling constants taken at saturation density />=0.16 fm^^ of the present work 
(first column) and from ref. ||2l| (second column, obtained from the Bonn A NN-potential) are 
compared to RMT values |l^ (third column). The meson masses used to derive the coupling 
constants from the results of Fig. 2 are displayed in the last column. 



18 



